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SUMMARY 


Forced  oscillation  of  a  two-dimensional  airfoil  with  attached  and 
separated  flows  is  investigated  using  nonlinear  unsteady  aerodynamics  for 
pitching  motion  derived  from  a  time  synthesization  technique  utilizing 
oscillatory  loop  data  determined  experimentally.  Both  one-  and  two-degree- 
of-freedom  oscillations  are  considered.  The  structural  dynamic  equations  of 
motion  are  integrated  by  a  time  marching  finite  difference  scheme.  The 
airfoil  response  is  examined  for  different  values  of  spring  stiffness  and 
magnitudes  of  externally  applied  moment.  For  two-degree-of-freedom 
vibration,  only  small  plunge  amplitude  is  considered  and  the  aerodynamic 
loads  are  approximated  by  the  superposition  of  nonlinear  terms  due  to 
pitch  and  linear  terms  due  tc  plunge.  The  presence  of  a  small  amplitude 
plunging  motion  increases  the  pitch  amplitude  slightly  for  attached  flow, 
while  a  decrease  in  pitch  amplitude  is  predicted  for  separated  flow. 


RESUME 


On  examine  l’oscillation  forcee  d’un  profil  de  voilure  bidimensionel 
dans  des  ecoulements  de  contact  et  separe  a  partir  de  donnees  d’aerodyna- 
rnique  non  lineaire  instable  sur  le  mouvement  de  tangage  produites  par  une 
technique  de  synthetisation  du  temps  basee  sur  une  boucle  oscillatoire 
etablie  experimentalement.  On  etudie  les  oscillations  a  un  et  a  deux  degres 
de  liberte.  Les  equations  dynamiques  structurales  du  mouvement  sont 
integrees  par  une  methode  d’avancement  du  temps  aux  differences  finies.  La 
reponse  du  profil  est  examinee  pour  differentes  valeurs  de  raideur  d’un 
ressort  et  de  moment  externe.  Pour  les  vibrations  a  deux  degres  de  liberte, 
seule  l’amplitude  des  faibles  plongeons  est  consideree  et  les  charges  aero- 
dynamiques  sont  approchees  par  la  superposition  de  termes  non  lineaires  de 
tangage  et  de  termes  lineaires  de  plongeon.  La  presence  d’un  mouvement  de 
plongeon  de  faible  amplitude  augmente  legerement  l’amplitude  du  tangage 
pour  l’ecoulement  de  contact,  tandis  qu’on  prevoit  une  diminution  de 
l’amplitude  du  tangage  pour  lecoulement  separe. 
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SYMBOLS 

Definition 
pitch  rate 

pitch  rate  at  dynamic  stall 

non-dimensional  distance  measured  from  airfoil  mid-chord  to  elastic  axis 
aerodynamic  sectional  static  pitching  moment  slope  at  zero  angle-of-attack 
aerodynamic  sectional  static  normal  force  curve  slope 
semi-chord  of  airfoil 

sectional  aerodynamic  pitching  moment  coefficient 

sectional  aerodynamic  static  pitching  moment  coefficient 

empirical  constants  in  Equation  (4) 

empirical  constants  in  Equation  (5) 

empirical  constants  in  Equation  (6) 

sectional  aerodynamic  normal  force  coefficient 

sectional  aerodynamic  static  force  coefficient 

chord 

natural  frequency  of  pitch  oscillation 
plunge  displacement 
reduced  frequency,  cub/U 
free  stream  Mach  number 
mass  of  airfoil  per  unit  span 
external  applied  force 
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Definition 

empirical  coefficients  in  Equation  7 
external  applied  moment 
amplitude  of  external  applied  moment 
empirical  constants  in  Equation  13 
Reynolds  number 

radius  of  gyration  about  elastic  axis 
time 

time  when  dynamic  stall  occurs 
free  stream  velocity 
defined  in  Equation  17 

non-dimensional  distance  measured  from  elastic  axis  to  centre  of  mass 
pitch  angle 

angle  of  attack  at  dynamic  stall 

mean  angle  of  attack 

angle  of  attack  at  reattachment 

angle  of  attack  at  static  stall 

angle  of  attack  when  vortex  leaves  trailing  edge 

defined  in  Equation  (1) 

value  of  at  dynamic  stall 

v'i'm- 

empirieal  constant  taken  as  0.18 
defined  in  Equation  (11) 
defined  in  Equation  (12) 
empirical  constant 


viscous  damping  for  pitching  motion 
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Symbol  Definition 

viscous  damping  for  plunging  motion 
ju  airfoil-air  mass  ratio  ,  =  m/7rpb" 

£  non-dimensional  displacement 

p  air  density 

r  non-dimensional  time 

Tm  non-dimensional  time  measured  from  instant  of  stall  onset 

rmt  defined  in  Equation  (6) 

0^  phase  angle  of  pitch  oscillation 

0C  defined  in  Equation  (3) 

to  frequency 

uncoupled  natural  frequency  of  plunging  motion 
coQ  uncoupled  natural  frequency  of  pitching  motion 

to  ratio  tJt/tJa 

Subscripts 


o 


initial  conditions 


FORCED  OSCILLATION  OF  A  TWO-DIMENSIONAL  AIRFOIL  WITH  NONLINEAR 

AERODYNAMIC  LOADS 


1.0  INTRODUCTION 

The  dynamic  response  of  a  two-dimensional  airfoil  to  external  oscillatory  forces  or  moments 
using  linear  aerodynamic  loads  derived  from  incompressible  flow  is  a  well  known  subject  (Refs.  1,2). 
For  more  complex  aerodynamics  as  in  transonic  flow,  numerical  time  marching  techniques  have  been 
used.  The  first  study  was  reported  by  Ballhaus  and  Goorjian  (Ref.  3)  who  carried  out  the  aeroelastic 
response  of  a  NACA64A006  airfoil  with  a  single-degree-of-freedom  in  pitch  at  transonic  speeds. 
Extensions  of  this  procedure  to  two-  and  three-degree-of-freedom  have  been  reported  by  Rizzetta 
(Ref.  4)  and  Yang  and  Chen  (Ref.  5).  Only  a  linear  response  was  treated  by  these  authors. 

In  aeroelastic  studies,  there  are  potentially  many  sources  of  nonlinearities  present,  but  the 
most  commonly  encountered  ones  are  those  having  structural  or  aerodynamic  origin.  Existing  tech¬ 
niques  of  analysing  dynamic  response  based  on  linear  vibration  theory  are  not  applicable.  Numerical 
methods  are  the  obvious  choice  in  solving  non-linear  vibration  problems  since  they  do  not  suffer  from 
the  limitations  of  perturbation  theory.  However,  they  have  received  little  attention  until  recently,  and 
this  study  shows  the  usefulness  of  the  numerical  approach  in  dealing  with  aerodynamic  nonlinearities 
associated  with  stalled  and  unstalled  airfoils.  This  problem  has  hitherto  not  been  amenable  to  theoret¬ 
ical  analysis.  A  prediction  method  of  the  dynamic  response  is  useful  in  providing  information  on  the 
sequence  of  events  occurring  on  the  airfoil  during  a  cycle  of  oscillation.  There  are  many  applications 
for  such  a  method,  for  example,  in  wind  tunnel  tests  of  oscillating  airfoils  at  high  incidence  or  large 
amplitude  oscillations.  In  helicopter  rotor  blades  design,  the  method  can  predict  unsteady  airloads 
and  deflections  of  the  blades  in  forward  flights  or  manoeuvring  operations. 

There  exists  a  number  of  numerical  time  marching  techniques  developed  for  finite  element 
linear  structural  analysis.  Among  the  most  commonly  used  time  integration  schemes  are  the  explicit 
central  difference  technique  and  Houbolt’s,  Wilson’s  and  Newmark’s  methods  (Refs.  6,  7).  Except  for 
the  explicit  scheme,  the  other  three  methods  are  unconditionally  stable.  Higher  order  schemes  are 
also  discussed  in  Reference  8  and  they  are  conditionally  stable.  Provided  that  care  is  taken  to  choose  a 
time  step  sufficiently  small  to  ensure  the  highest  mode  considered  in  a  vibrating  system  does  not 
diverge,  higher  order  schemes  are  more'  efficient  in  terms  of  computation  time.  However,  for  a  system 
with  few  structural  components  or  the  number  of  vibration  modes  is  small,  higher  order  methods 
such  as  the  eighth  order  scheme  reported  in  Reference  8  does  not  offer  any  distinct  advantage  over 
Houbolt's  scheme  which  is  simpler  and  less  cumbersome  to  use. 

The  use  of  numerical  time  integration  techniques  to  study  nonlinear  vibration  in  one-degree- 
of-freedom  was  reported  in  Reference  8.  The  nonlinearities  considered  was  that  of  a  cubic  spr  m.  The 
numerical  results  agree  very  well  with  analytical  predictions  derived  from  perturbation  theory,  and 
in  addition,  give  more  information  on  the  behaviour  of  the  system  to  initial  conditions  which  the 
analytical  method  fails  to  provide. 

In  this  study.  Houbolt's  (Ref.  7)  scheme  is  used  to  investigate  the  forced  oscillation  of  an 
airfoil  in  one-  and  (wo-ilegree-of-freedom.  For  the  one-degree-of-freedom  motion  case,  only  pitch 
oscillation  is  investigated  and  both  stalled  and  unstalled  flow  are  considered.  There  are  a  number  of 
studies  l  Refs.  9-14  )  on  methods  of  predicting  dynamic  stall  and  unsteady  airloads  on  two-dimensional 
airfoils  with  harmonic  pitching  motion.  The  most  sophisticated  one  is  given  by  Bielawa  et  al.  (Ref.  14  t 
using  experimental  oscillatory  loop  data  to  generate  synthesized  data  in  the  time  domain.  This  method 
is  used  m  (he  present  study  since  it  conveniently  generates  the  aerodynamic  loads  at  each  time  step 
m  the  integration  of  the  structural  dynamic  equations  of  motion. 

In  the  t  wo-degree-of-freedom  vibration  of  the  airfoil,  the  plunge  amplitude  is  assumed  to 
be  small  Hie  aerodynamic  loads  are  then  approximated  by  the  superposition  of  nonlinear  terms  due 
to  pitch  and  linear  terms  due  to  plunge,  since  empirical  relations  for  both  large  amplitude  pitch  and 


plunge  motions  are  not  available.  A  comparison  with  one-degree-of-freedom  response  results  gives  the 
effect  of  a  small  plunging  motion  for  cases  with  attached  or  separated  flow  on  the  airfoil. 


2.0  ANALYSIS 

2.1  Empirical  Representation  of  LInsteady  Aerodynamic  Loads  of  Stalled  and  Unstalled  Airfoils 

In  Reference  14  an  empirical  method  to  determine  the  aerodynamic  loads  is  given  in  the 
time-domain  using  data  obtained  from  oscillating  airfoil  experiments.  The  expressions  for  the  normal 
force  and  pitching  moment  coefficients  are  quite  general  and  valid  for  both  stalled  and  unstalled 
airfoils.  For  a  given  airfoil  shape,  Mach  number  and  Reynolds  number,  the  dynamic  characteristics  of 
the  airfoil  depend  on  the  mean  angle  of  attack,  the  frequency  and  amplitude  of  oscillations.  In  the 
case  of  dynamic  stall,  it  is  assumed  that  a  vortex  develops  near  the  leading  edge  when  the  static  stall 
angle  is  exceeded.  As  the  angle  of  incidence  a  increases,  the  vortex  detaches  from  the  leading  edge 
and  convects  downstream  near  the  surface  until  it  leaves  the  trailing  edge.  The  airfoil  remains  stalled 
until  o  drops  significantly  for  reattachment  of  flow  to  occur. 

Bielawa  et  al.  (Ref.  14)  define  a  parameter  aw  which  accounts  for  the  time  history  effects 
of  the  change  in  a  and  is  given  by 
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where  r  is  the  non-dimensional  time 
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=  V  1  -M-,  a (0)  is  the  initial  angle  of  attack  at  time  r  =  0,  and  <£c(r,M)  is  the  compressibility 
corrected  Wagner  function  which  is  written  as 


0c(r, M)  =  -  1  -  0.1(i5e'"  "455  r')"  -  0.335e  0  3  rJ" 


(3) 


The  equations  in  the  rest  of  this  section  are  essentially  the  same  as  those  obtained  from 
Reference  14.  They  are  included  here  since  they  are  used  in  the  numerical  finite'  difference  scheme  to 
be  described  in  the  next  section.  For  a  detailed  description  of  the  time  synthesization  technique,  the 
derivation  given  by  Bielawa  (Ref.  14)  should  be  referred  to. 

The*  dynamic  stall  and  reattachment  angles  in  terms  of  the  static  stall  angle  «ss.  the  pitch 
rate  A  and  are  defined  as 
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where  c.  CAln.  CWnl.  C  AK  and  R  are  empirical  constants  to  be  determined.  The  following  relationship 
is  used  to  predict  the  time  it  takes  for  the  vortex  to  travel  from  the  leading  edge  to  the  trailing  edge  of 
the  airfoil : 


(12) 


The  pitching  coefficient  is  expressed  as 


cm  =  Cms  (<*-Ao2)  +  at,M Ao2  +  Q,A  +  Q2aw  +  Q3 


+  Q4  |aw|  +  Q5 ®  1  +  ^6 Ao2  +  Q7«DmADmrm  (13) 

where  a.,,  is  the  static  pitching  moment  slope  at  zero  angle  of  attack.  For  unstalled  airfoils,  the  last 
UM 

three  terms  in  Equation  (13)  are  zero. 

2.2  Two-Degree-of-Freedom  Motion  of  a  2-D  Airfoil 

Figure  1  shows  the  notations  used  in  the  analysis  of  a  two-degree-of-freedom  motion  of  an 
airfoil  oscillating  in  pitch  and  in  plunge.  The  plunging  deflection  is  denoted  by  h,  positive  in  the  down 
ward  direction,  a  is  the  pitch  angle  about  the  elastic  axis,  positive  with  the  nose  up.  The  elastic  axis 
is  located  at  a  distance  altb  from  the  midchord,  while  the  mass  centre  is  located  at  a  distance  xQb 
from  the  elastic  axis.  Both  distances  are  positive  when  measured  towards  the  trailing  edge  of  the 
airfoil.  The  aeroelastic  equations  of  motion  have  been  derived  by  Fung  (Ref.  2)  and  can  be  written  as 
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where  $  =  h/h  is  the  non-dimensional  displacement,  m  is  the  mass  per  unit  span  of  the  airfoil  and 


cjf,  tJ0  are  the  uncoupled  plunging  and  pitching  natural  frequency  respectively,  f.  and  are  the 
viscous  damping  ratio  for  plunge  and  pitch  motion  respectively,  ia  is  the  radius  of  gyration  about  the 
elastic  axis,  CN  and  CM  are  the  normal  force  and  pitching  moment  coefficients,  and  P(r)  and  Q(t)  are 
the  external  applied  force  and  moment  respectively.  U*  is  defined  as 

U 

u*  =  - -  (17) 


For  given  P(t)  and/or  Q(r),  Equations  (14)  and  (15)  can  be  solved  for  the  forced  oscillation 
of  the  airfoil  in  pitch  and  plunge  provided  that  CN(r)  and  CM(r)  are  known.  In  linear  analysis  for 
small  oscillations  of  the  airfoil,  superposition  of  CN  and  CM  for  pitch  and  plunge  motions  is  permis¬ 
sible  in  determining  the  total  force  and  moment  coefficients.  However,  when  the  aerodynamic  loads 
are  nonlinear,  these  coefficients  have  to  be  determined  for  combined  motions  in  pitch  and  plunge. 

The  two-degree-of-freedom  oscillating  airfoil  studied  in  this  report  has  large  stiffness  in  the  plunge 
motion.  In  the  limit  when  co  is  infinite,  the  motion  degenerates  to  that  of  a  one-degree-of-freedom 
system.  Since  the  plunging  motion  is  assumed  to  be  small,  the  values  of  CN  and  CM  are  predominantly 
due  to  pitching  motion.  The  contributions  due  to  plunge  are  added  on  from  linear  aerodynamic  theory 
and  hence  they  should  be  treated  as  corrections  which  are  only  approximations  to  an  otherwise 
extremely  complex  situation. 

2.3  Finite  Difference  Scheme 

Houbolt’s  (Ref.  7)  implicit  method  is  used  in  the  present  analysis  even  though  more 
accurate  higher  order  schemes  are  available  (Ref.  8).  In  this  case,  the  derivatives  at  time  r  +  A r  are 
replaced  with  backward  difference  formulas  using  values  at  three  previous  points.  For  example, 

“nar  =  “~7  (2ar+Ar  -  +  4«r-AT  ~  ar -2At)  +  (At2)  (18> 
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and 

® r  +  A T  =  TT“  (llaT+Ar  "  18ftr  +  9°7-At  "  2°V-2At)  +  0(Ar3)  (19) 


Similar  expressions  can  be  written  for  £r  +  Ar  and  £T+Ar.  In  difference  form.  Equations  (14)  and  (15) 
can  be  written  as 

EW  +  Vor+Ar  =  T  (20) 

M£T+Ar  +  Iar+Ar  =  U  (21) 

The  coefficients  E,  V,  M,  1  and  the  terms  T  and  U  on  the  RHS  of  Equations  (20)  and  (21)  are  given 
in  the  Appendix. 

2.4  Starting  Procedure 

Houbolt’s  scheme  requires  values  of  a  and  £  at  times  r-2Ar,  r-Ar  and  r  in  order  to 
determine  their  values  at  r+Ar.  At  time  r  =  0  a  special  starting  procedure  is  required.  Writing 
Equations  (14)  and  (15)  at  r  =  0  and  solving  for  q0  and  £„  gives 
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where  p0  and  r0  can  be  obtained  from  the  Appendix,  and  the  initial  conditions  a0,  a0,  £0  and  are 
known.  A  Taylor  series  is  then  used  to  obtain  the  following 


At2  .. 


c*_Ar  =  ao  ~  Ara0  +  —  o0  +  0(Ar3) 


a0  +  Ard0  +  —  a0  +  0(Ar3) 


with  similar  expressions  for  £_Ar  and  fcAr.  For  the  next  step,  Houbolt’s  scheme  can  be  used  since 
a-Ar>  ao’  aAr>  £-at>  So  anc*  Sat  are  known.  The  accuracy  of  the  numerical  method  is  0(Ar4)  on  each 
step  while  Equations  (24)  and  (25)  limit  the  accuracy  to  0(At3).  A  starting  accuracy  higher  than 

,  2n  , 

0(At  )  is  not  necessary  since  the  error  per  cycle  in  the  numerical  scheme  is  —  0(Ar  )  (Ref.  8). 

cu 

3.0  RESULTS  AND  DISCUSSIONS 

3.1  Synthesized  Data  for  a  Vertol  Modified  NACA0012  Airfoil 

In  Reference  15  two-dimensional  oscillatory  airfoil  test  data  sets  for  pitching  motion  are 
given  for  a  Vertol  Modified  NACA0012  airfoil.  In  this  report  the  synthesized  data  are  only  shown 
for  M  =  0.6  and  Rc  =  6.2  X  106.  Equations  (4)  to  (6)  predict  the  stall  events  and  the  coefficients  in 
these  equations  are  determined  empirically.  The  force  and  moment  coefficients  are  obtained  from 
Equations  (7)  and  (13)  by  a  curve  fitting  procedure  using  data  loops  for  both  unstalled  and  stalled 
conditions.  In  this  particular  example,  13  data  loops  are  used  for  CN  and  14  data  loops  for  CM.  The 
coefficients  in  the  two  equations  are  obtained  by  a  minimization  procedure  given  by  Powell  (Ref.  16) 
Usually  in  each  cycle  of  oscillation  of  the  airfoil,  the  loop  is  divided  into  600  time  steps.  The  steps 
are  adjusted  so  that  the  spacings  are  reduced  in  regions  where  large  changes  in  aerodynamic  charac¬ 
teristics  occur. 

In  Table  1,  the  empirical  coefficients  in  Equations  (4)  to  (6)  and  those  for  C’N  and  CM  at 
M  =  0.62  and  Re  =  6.2  X  106  are  given.  The  comparisons  between  the  synthesized  and  test  data  are 
shown  in  Figures  2  and  3.  The  correlation  is  very  simila'  to  that  obtained  by  Bielawa  et  al.  (Ref.  11 ) 
which  is  considered  to  be  good  compared  to  other  empirical  formulations. 


3.2  Forced  Oscillation  for  Pitching  Motions 


For  very  large  values  of  cu,  Equations  (14)  and  (15)  reduce  to  that  for  pitch  oscillation.  The 
forcing  function  in  Equation  (15),  that  is,  Q(T)/mU3  can  be  written  as  Q()  sin  kr  where  the  reduced 
frequency  k  =  cob/U. 

Figure  4a  shows  the  angular  displacement  from  the  mean,  CN  and  CM  for  the  first  five 
cycles  after  an  external  moment  has  been  applied  at  r  =  0.  The  pitch  axis  is  at  the  1/4  chord,  and  the 
mean  angle  of  attack  am  is  0.2°.  The  initial  displacement  and  angular  velocity  of  the  airfoil  are  zero. 
The  driving  frequency  corresponds  to  that  for  a  value  of  k  =  0.165.  The  amplitude  of  the  applied 
moment  is  Qu  =  0.8  X  10'3,  f a  =  0  and  co/oja  =  0.9.  The  airfoil  has  the  following  properties:  p  =  100, 
ra  =  0.5,  xa  -  0.25  and  a,,  =  -0.5.  The  flow  around  the  airfoil  is  attached  at  all  times  and  the 
oscillations  reach  a  steady  states  in  four  or  five  cycles.  The  results  between  15  to  20  cycles  are  shown 
in  Figure  4b  and  it  is  seen  that  they  are  practically  the  same  as  those  at  the  fifth  cycle.  In  all  compu¬ 
tations,  a  time  step  equal  to  1/128  of  a  cycle  is  used  and  found  to  give  sufficiently  accurate  results. 

At  larger  mean  angle  of  attack  when  the  airfoil  stalls,  it  usually  takes  a  few  more  cycles 
before  a  steady  state  is  reached.  Figures  5a  and  5b  show  the  first  five  and  the  15  to  20  cycles  for 
am  -  7.48°,  Qt)  =  0.5  X  10~3,  cj/<ua  =  1.4  and  k  =  0.165. 

Figures  6  and  7  show  the  amplitude  and  phase  of  an  unstalled  airfoil  at  otm  =  0.2°,  M  =  0.6 
for  three  values  of  Q0.  The  driving  frequency  is  kept  constant  at  a  value  of  k  =  0.165,  and  varying 
to/coQ  is  equivalent  to  changing  the  stiffness  of  the  torsional  spring  constant.  Since  f  =  0,  the 
damping  is  solely  from  the  aerodynamics.  The  phase  curves  for  Q0  =  0.5  and  0.8  X  10_3are  shifted 
upward  by  20°.  These  two  figures  are  very  similar  to  those  for  one-degree-of-freedom  system  with 
viscous  damping. 

With  the  same  value  of  k  but  increasing  am  to  7.48°,  the  amplitude  and  phase  curves 
are  shown  in  Figures  8  and  9  for  six  values  of  the  amplitude  of  the  external  driving  moment  Qu. 

Again,  the  phase  curves  are  shifted  upwards  by  20°.  Except  for  curve  ‘6’  with  the  smallest  value  of 
Q(,  =  0.1  X  10'3,  the  other  five  cases  exhibit  breaks  in  the  amplitude  and  phase  versus  cj/cua  curves. 
Starting  with  small  values  of  o;/coa,  the  flow  over  the  airfoil  is  attached  until  a  value  of  co/u>a  is 
reached  where  after  many  cycles  of  computation  a  steady  condition  does  not  appear  to  exist.  The 
flow  changes  from  attached  to  separated  and  back  to  attached  and  back  and  forth  without  any 
definite  pattern.  Further  increase  in  ou/coa  will  result  in  a  steady  condition  with  separated  flow  over 
the  airfoil.  A  maximum  in  the  amplitude  of  oscillation  of  the  airfoil  is  reached  in  the  vicinity  of 
oj/wu  =  1  and  the  amplitude  and  phase  curves  behave  like  those  for  a  linear  oscillator.  However, 
upon  increasing  co/coa,  breaks  in  the  curves  are  again  detected  and  for  the  two  smaller  values  of  Q0, 
i.e.  Q0  =  0.25  and  0.5  X  10'3,  there  is  a  small  region  of  unsteadiness  where  the  flow  does  not  settle 
either  to  the  attached  or  separated  conditions,  but  beyond  which  the  curves  reach  a  steady  condition 
again  with  the  airfoil  oscillating  in  the  unstalled  state.  For  the  larger  values  of  Qu,  no  steady  conditions 
can  be  reached.  The  failure  to  reach  steady  oscillations  when  the  flow  changes  from  attached  to 
separated  or  vice  versa  in  those  regions  where  the  breaks  occur  is  probably  due  to  the  method  of 
calculating  CN  and  CM  from  Equations  (7)  and  (13).  At  each  step  in  the  numerical  finite  difference 
scheme,  the  local  values  of  a  and  d  are  used  to  compute  the  local  pitch  rate  A  and  aw.  These  are 
then  substituted  into  Equation  (4)  to  evaluate  a  value  of  a1)m.  Depending  on  whether  a  is  greater  or 
less  than  aj)in,  the  flow  is  taken  to  be  either  separated  or  attached  accordingly.  In  the  first  transition 
region,  the  amplitude  a  is  initially  smaller  than  a()||1  and  the  flow  is  attached.  As  the  amplitude  grows, 
«  will  exceed  a1)ni  and  the  flow. separates.  Because  of  the  ensuing  increase  in  damping  a  then  decreases 
and  the  flow  becomes  attached  again  with  a  smaller  value  of  damping.  The  value  of  a  then  starts  to 
increase  and  the  cycle  repeats  itself.  The  same  phenomenon  also  occurs  in  the  second  transition  region. 
To  eliminate  this  oscillation  between  attached  and  separated  flows,  a  different  scheme  to  fix  the  state 
of  the  flow  at  transition  has  to  be  devised. 

Instead  of  holding  k  constant,  Figures  10  and  11  show  the  amplitude  and  phase  by  varying 
the  frequency  of  the  external  moment  for  the  airfoil  at  cvni  =  7.48°,  Q„  =  0.5  X  10"3  and  natural 
frequencies  of  24,  48,  64  and  80  Hz.  These  curves  show  the  same  characteristics  as  those  in  Figures  8 


and  9.  The  regions  of  unsteadiness  where  the  flow  changes  from  either  attached  to  separated  flows  or 
vice  versa  are  much  larger  for  low  natural  frequencies. 

The  results  given  so  far  are  for  =  0.  For  small  values  of  the  damping  ratio  fa,  Figure  12 
shows  a  typical  modification  of  the  response  curve  at  k  =  0.165  and  Q()  =  0.8  X  10'3.  Viscous 
damping  tends  to  decrease  the  amplitude  and  increase  slightly  the  region  of  unsteadiness  where  the 
flow  oscillates  between  attached  and  separated  states. 

3.3  Forced  Oscillations  for  Pitching  and  Plunging  Motions 

The  empirical  relations  given  in  Section  2.1  for  the  unsteady  aerodynamic  loads  are  for 
pitching  motion  only.  When  the  airfoil  is  oscillating  in  two-degree-of-freedom  with  large  amplitudes 
of  motion,  a  suitable  method  of  representing  the  nonlinear  aerodynamics  has  to  be  formulated. 
Experimental  oscillatory  data  for  the  synethesization  method  have  to  include  cases  for  various  com¬ 
binations  of  pitching  and  plunging  amplitudes.  These  data  are  difficult  to  obtain  and  are  not  presently 
available. 

The  present  investigation  considers  the  case  of  a  two-degree-of-freedom  motion  with  large 
amplitude  in  pitch  but  small  amplitude  in  plunge.  In  other  words,  the  stiffness  of  the  spring  for 
plunging  motion  is  kept  large  so  that  to  in  Equations  (14)  and  (15)  is  large.  The  aerodynamic  loads 
are  then  given  by  the  sum  of  two  terms:  pitching  motion  from  Equations  (7)  and  (13)  and  plunging 
motion  from  linear  aerodynamics  given  in  Reference  1  using  the  indicial  lift  and  moment  functions 
at  M  =  0.6  determined  by  Mazelsky  and  Drischler  (Ref.  17).  This  formulation  of  the  aerodynamic 
loads  is  only  approximate  but  it  is  used  in  this  study  to  give  some  idea  of  the  effect  of  a  plunge  degree 
of  freedom  motion  with  small  amplitude  on  the  pitching  motion  of  the  airfoil  driven  by  an  externally 
applied  moment. 

For  an  unstalled  airfoil  with  am  =  0.2°,  fa  =  =  0,  Figure  13  shows  the  pitch  and  plunge 

motions  and  the  aerodynamic  coefficients  when  steady  conditions  are  reached.  The  value  of  the 
torsional  natural  frequency  is  fa  =  64  Hz,  co  =  2,  Qt,  =  0.5  X  10'3  and  the  results  are  given  for 
co/toa  =  2.2  for  cycles  35  to  40.  In  all  computations,  a(0)  =  a(0)  =  £(0)  =  £(0)  =  0.  The  effects  of 
fa  and  to  on  the  amplitude  response  with  variation  in  co/coQ  are  shown  in  Figures  14  to  19.  The 
amplitude  of  the  pitching  motion  increases  with  decreasing  to  for  the  three  values  of  fQ  considered, 
that  is,  for  the  same  value  of  the  applied  moment  (Q0  =  0.5  X  10~3),  the  presence  of  a  plunging  motion 
increases  the  pitch  amplitude  slightly.  The  second  peak  in  the  vicinity  of  to  is  usually  quite  small.  For 
the  £  response  curves  shown  in  Figures  15,  17  and  19,  the  second  peak  is  comparable  to  and  in  some 
cases  larger  than  the  first.  However,  it  is  not  strongly  dependent  on  fo  and  its  magnitude  changes 
only  slightly  with  increasing  f  which  is  quite  unlike  the  first  peak  where  the  amplitude  drops  very 
rapidly  as  the  torsional  natural  frequency  is  increased. 

When  the  airfoil's  mean  angle  of  attack  is  increased  to  cvin  =  7.48°,  the  response  is  similar 
to  that  for  one-degree-of-freedom  motion  discussed  in  the  previous  section  (see  Fig.  10).  The  effect 
of  a  plunge  degree-of-freedom  on  the  amplitude  of  the  pitching  motion  is  shown  in  Figure  20  for 
f0  =  64  Hz,  Q„  =  0.5  X  10‘3  and  to  =  10,  4  and  3  respectively.  For  small  values  of  co/to,,  the  flow  is 
attached  and  the  transition  from  attached  to  separated  flow  occurs  around  to/toa  =  1.  to  does  not 
appear  to  have  much  effect  on  the  range  of  the  unsteady  region  where  the  airfoil  oscillates  between 
attached  and  separated  flow  conditions.  Unlike  Figures  14,  16  and  18  which  show  a  slight  increase 
in  the  amplitude  of  o  as  to  decreases.  Figure  20  indicates  the  opposite,  that  is,  for  a  stalled  airfoil, 
introducing  a  plunge  degree-of-freedom  will  decrease  the  amplitude  of  the  pitch  motion  for  the  same 
magnitude  of  the  externally  applied  moment .  As  to  to  ,  increases,  the  flow  becomes  attached  again. 
For  large  values  of  to  (to  =  10  in  this  particular  case),  a  steady  condition  is  reached  after  approxi¬ 
mately  twenty  cycles  of  forced  oscillation.  However,  decreasing  to  results  in  failure  to  achieve  a 
steady  condition  even  over  50  cycles  of  computation.  Tin-  flow  remains  attached  but  the  amplitudes 
of  the  pitch  and  plunge  motions  are  scattered  as  shown  m  curves  2  and  3  of  Figure  20.  There  is  no 
plausible  explanation  for  this  anomaly  at  present 


4.0  CONCLUSIONS 


The  empirical  relations  used  to  represent  nonlinear  aerodynamic  force  and  moment 
coefficient  of  both  stalled  and  unstalled  airfoils  give  fairly  good  results  compared  with  the  original 
loop  data  from  which  they  are  derived.  They  are  formulated  for  pitching  motion  only  and  modifi¬ 
cations  have  to  be  made  if  they  are  to  be  used  for  plunging  or  combination  of  plunging  and  pitching 
motions.  The  synthesized  results  can  be  used  to  derive  data  for  any  initial  angle  of  attack  of  the 
airfoil,  pitch  amplitude  and  frequency  of  oscillation,  but  in  the  present  formulation  are  restricted  to 
the  same  airfoil  and  Mach  number  from  which  the  experimental  oscillatory  aerodynamic  data  are 
obtained. 


In  one-degree-of-freedom  forced  oscillation  in  pitch,  steady  conditions  for  unstalled  flow 
can  be  reached  after  a  few  cycles  of  motion  of  the  airfoil  starting  from  rest  initially.  For  a  stalled 
airfoil,  usually  a  few  more  cycles  of  computations  are  required  to  reach  steady  state.  In  those  cases 
where  the  flow  changes  from  attached  to  separated  or  vice  versa  the  transition  regions  are  very 
unsteady  and  the  flow  is  unsettled  between  the  stalled  and  unstalled  conditions  without  reaching  a 
steady  state,  irrespective  of  the  number  of  computation  cycles  carried  out.  Neither  does  a  decrease  in 
the  time  step  used  in  the  numerical  finite-difference  scheme  show  any  improvement.  This  is  probably 
due  to  an  inadequate  method  used  to  compute  the  dynamic  stall  angle  at  those  critical  forcing 
frequencies.  The  inclusion  of  a  viscous  damping  term  in  the  dynamic  equation  of  motion  does  not 
introduce  any  drastic  changes  in  the  behaviour  of  the  response  aside  from  a  decrease  in  amplitude 
and  slight  increase  in  the  unsteady  region  where  the  flow  oscillates  between  attached  and  separated 
states. 


For  two-degree-of-freedom  motion,  only  small  plunge  amplitude  is  considered  because  of 
the  approximations  used  in  computing  the  aerodynamic  loads.  The  presence  of  a  plunging  motion 
increases  the  pitch  amplitude  slightly  for  attached  flow,  while  a  decrease  in  the  pitch  amplitude  is 
predicted  for  separated  flow. 
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FIG.  5b:  ANGULAR  DISPLACEMENT  FROM  THE  MEAN,  CN  AND  CM 
FOR  THE  15-20  CYCLES  OF  FORCED  PITCH  OSCILLATION, 
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APPENDIX 


SUMMARY 


Forced  oscillation  of  a  two-dimensional  airfoil  with  attached  and 
separated  flows  is  investigated  using  nonlinear  unsteady  aerodynamics  for 
pitching  motion  derived  from  a  time  synthesization  technique  utilizing 
oscillatory  loop  data  determined  experimentally.  Both  one-  and  two-degree- 
of-freedom  oscillations  are  considered.  The  structural  dynamic  equations  of 
motion  are  integrated  by  a  time  marching  finite  difference  scheme  The 
airfoil  response  is  examined  for  different  values  of  spring  stiffness  and 
magnitudes  of  externally  applied  moment.  F'or  two-degree-of-freedom 
vibration,  only  small  plunge  amplitude  is  considered  and  the  aerodynamic 
loads  are  approximated  by  the  superposition  of  nonlinear  terms  due  to 
pitch  and  linear  terms  due  to  plunge.  The  presence  of  a  small  amplitude 
plunging  motion  increases  the  pitch  amplitude  slightly  for  attached  flow, 
while  a  decrease  in  pitch  amplitude  is  predicted  for  separated  flow. 


RESUME 


On  examine  l’oscillation  forcee  d’un  profil  de  voilure  bidimensionel 
dans  des  ecoulements  de  contact  et  separe  a  partir  de  donnees  d’aerodyna- 
mique  non  lineaire  instable  sur  le  mouvement  de  tangage  produites  par  une 
technique  de  synthetisation  du  temps  basee  sur  une  boucle  oscillatoire 
etablie  experimentalement.  On  etudie  les  oscillations  a  un  et  a  deux  degres 
de  liberte.  Les  equations  dynamiques  structurales  du  mouvement  sont 
integrees  par  une  methode  d’avancement  du  temps  aux  differences  finies.  La 
reponse  du  profil  est  examinee  pour  differentes  valeurs  de  raideur  d’un 
ressort  et  de  moment  externe.  Pour  les  vibrations  a  deux  degres  de  liberte. 
seule  l’amplitude  des  faibles  plongeons  est  consideree  et  les  charges  aero- 
dynamiques  sont  approchees  par  la  superposition  de  termes  non  lineaires  de 
tangage  et  de  termes  lineaires  de  plongeon.  La  presence  d'un  mouvement  de 
plongeon  de  faible  amplitude  augmente  legerement  l’amplitude  du  tangage 
pour  I’ecoulement  de  contact,  tandis  qu’on  prevoit  une  diminution  de 
l’amplitude  du  tangage  pour  l’ecoulement  separe. 


